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ON THE SHAPE OF CORE OVERSHOOTING IN STELLAR 
MODEL COMPUTATIONS, AND ASTEROSEISMIC TESTS 

Ehsan Moravveji * 1 


Abstract. Slowly pulsating B stars (SPB) and 7 Dor stars pulsate in 
high-order gravity (g-) modes. The frequencies of g-modes are sensitive 
to the detailed structure and evolution history of stars having convec¬ 
tive cores. Receding convective cores in OB-type stars leave behind a 
chemically inhomogenous > 0 radiative zone. Once a g-mode has 
radial nodes near the boundaries of these layers, the mode gets trapped 
and its period deviates from asymptotic period spacing. Careful study 
of such trapped modes allows constraining the extent of such layers 
by fitting individual pulsation frequencies. We employ 19 consecuitve 
dipole g-modes of a very rich Kepler SPB pulsator, KIC 10526294, to 
demonstrate the power of mode trapping in B-stars in studying the 
thermal and chemical stratification in the overshooting layer. 


1 Introduction 


Once OB-type stars reach the zero-age main sequence (ZAMS), they immediately 
develop fully mixed convective cores, which already attains the largest possible 
mass and extent at ZAMS. During the evolution, the drop in center hydrogen (X c ) 
is accompanied by the decrease in radiative temperature gradient V ra d- Conse¬ 
quently, the convective core shrinks monotonically. This leaves behind CNO en¬ 
riched material and a steep increase in hydrogen abundance outside the boundary. 
Therefore, = (d In fi/d In p) > 0; here, fi is the mean molecular weight, and 
p is the total pressure. The boundary layer mixing induced by the overshoot of 
convective plumes into the stable radiative layers can partially mix the hydrogen, 
helium and metals outside the core, and reshape the profile there. 

The problem of the boundary layer mixing is still not understood well. Re¬ 


cently, Viallet et al. (2015) proposed that the overshooting mixing depends on 
the local Peclet number - defined as the local ratio of the thermal to viscous 
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timescales - in the overshooting layer: Closer to the core boundary, the convective 
plumes are predicted to be quasi-adiabatic, increasing the size of the fully mixed 
core. Furhter away, the Peclet number drops steeply, and the overshooting layer 
becomes photon-dominated, and the mixing proceeds diffusively. These two in¬ 
fluence the shape of just ouside the core, and it propagates directly into the 
pulsation equations through the definition of the Brunt-Vaisala frequency 


^iv = #(v ad-V + ^V, 

tin 


( 1 . 1 ) 


where H p is the local pressure scale height, V = (dlnT/dlnp) is the actual tem¬ 
perature gradient V a d is the adiabatic temperature gradient, S = (d In p/d In T)p jIJL 
and (p = (d\np/d\np) p^t- The local thermal and compositional stratifications are 
encapsulated in V and in Eq. (|1.1|) . 


The Brunt-Vaisala frequency is explicitly present in the adiabatic and non- 

1989; Aerts et al 


adiabatic linearized oscillation equations (Unno et al. 


2010 ) 


Thus, even a slight modification to the by mixing (through overshooting or 
extra diffusive mixing) in the radiative part of the star influences the eigenfunc¬ 
tions and eigenfrequencies, and exhibits measurable fingerprints in period spacing 
A P = P n+ 1 — P n , where P n is the period of a dipole g-mode of radial order n. 
Asteroseismology of heat-driven g-mode pulsators opens a direct window to study 
the physical properties of the fully-mixed core, and the partially homogeneous 


layer on top of the core (Miglio et al. 2008). 


Recently, [Papics et al. (2014) studied the SPB star KIC 10526294, and iden¬ 


tified 19 dipole g-modes from the triplet structure around each mode. Mor avveji| 
et al. ( 2015| MAP 15) carried out a detailed forward seismic modelling of the dipole 
zonal frequencies, and |Triana et al. (2015) inferred the internal rotation profile by 
inverting its rotational splittings. We demonstrate how trapping of high radial or¬ 
der g-modes in the V^-layer is used to confine the size of this layer. We also argue 
that our current treatment of the overshoot mixing requires a major revision. 


2 Simplified picture of mixing in B stars 


MAP15 compared the exponentially decaying ( D ov = D conv exp[—2(r— r core ) / f ov H p )]) 
and the step-function ( D ov = D conv for r core <r< a ov H p ) overshooting prescrip¬ 
tions in their extensive seismic modelling of KIC 10526294. Note that a ov ~ 10/ ov . 


2013 


2000 


The mixing induced by convective overshoot is assumed difhusive (Zhang 
|Viallet et al. 2015) with an exponentially decaying behaviour (Herwig 
(Figj^t) or step-function behaviour (Fig[]J>). We also assume that the thermal 
stratification in the overshoot layer is purely radiative V = V ra d- Clearly, N^ v 
(red dashed) shows a steeper rise in the step-overshoot model. Fig. lc shows the 
observed versus modelled AP, where the first five modes agree very well with ob¬ 
servations, and the other thirteen deviate significantly from the fit. We argue the 
reason for this is our reasonable ability in confining the width of the //-gradient 
layer, and at the same time our poor understanding of the thermal and chemical 
stratification in the overshooting layer. 


















































The overshooting properties of KIC 10526294 
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Fig. 1 . The convective (blue), overshooting (grey) and extra diffusive mixing (orange) 
profiles in the best models of KIC 10526294, using exponential (a) and step-function (b) 
overshooting prescriptions. See Table 4 in MAP 15 for the mass and extent of each zone. 
The dashed lines show the Brunt-Vaisala N^ v (Eq |l.l| ) profiles, (c) The observed (grey) 
versus modelled (black) period spacing from the best seismic model (with / ov = 0.017). 


3 Mode trapping 


Fig. [2 presents the behaviour of the rotation kernels K n ^ (defined in Aerts et al. 
(2010), Eq. 3.356) versus enclosed mass. The bottom panels show the run of the 
Brunt-Vaisala frequency Nbv, the Lamb frequency Si and the mode frequency 
in logarithmic scale. The /i-gradeint zone corresponds to the local peak in Nbv 
around 0.85 M 0 . The shortest period mode (left) is perfectly trapped in the /i- 


gradient layer, and has a sizable amplitude there (Dziembowski and Pamyatnykh 


1991). The longest period mode (right) still has a large amplitude in the /i-gradient 
layer, in addition to a significant amplitude in the diffusive oveshooting layer below 
the //-gradient zone. Thus, the eigenfunctions of the latter mode are sensitive to the 
treatment of overshooting (FigsjTJi andjlj)). The frequencies and period spacings 
assosicated to highest-order g-modes of KIC 10526294 (Fig.JTj^) do not perfectly 
match the observation. This manifests that our underlying assumptions about the 
thermal and compositional structure of the overshooting layer build in ID stellar 
structure codes need a major revision. 

The current assumption in ID stellar models is that the temperature gradient 
V changes abruptly from fully adiabatic in the convective core to fully radiative 
just above the convective core, where the convective plumes overshoot into the 
radiative regions. However, this picture is not supported by 3D simulations of 
Viallet et al. (2013| (although performed for the conditions of deep convection in 
red giant envelopes). It is possible that the behaviour of the temperature gradient 
in the overshooting layer is quasi-adiabatic, up to some pressure scale heights from 
the core boundary, where it efficiently mixes the materials and adds to the size of 
the fully mixed core. Such a prescription is similar to the classical step-function 
overshoot. An extension of this quasi-adiabatic layer can be photon-dominated, 
and an exponentially decaying diffusive mixing can be applicable there. There¬ 
fore, a more complicated ID stellar evolution models should be built, hopefully 
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Fig. 2. The rotation kernels for the shortest period mode n — —14 (Left) and for the 
longest period mode n — —32 (Right) chosen from the best model of KIC 10526294. The 
profile of the Brunt-Vaisala frequency IVbv and the Lamb frequency are shown to guide 
the eyes. The convective core is highlighted in blue. 


with improved treatment of convection (Zhang, 2013 Arnett et ah, 2015) to pro¬ 


vide better inputs for asteroseismic inferences. Indeed, the one-to-one frequency 
matching with observations can discriminate between the most plausible convec¬ 
tive/overshoot prescription. 
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